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ABSTRACT 

We perform deep 1.8 cm radio continuum imaging towards thirteen protostellar regions se- 
lected from the Spitzer c2d small clouds and cores programme at high resolution (25") in 
order to detect and quantify the cm-wave emission from deeply embedded young protostars. 
Within these regions we detect fifteen compact radio sources which we identify as radio pro- 
tostars including two probable new detections. The sample is in general of low bolometric 
luminosity and contains several of the newly detected VeLLO sources. We determine the 
1.8 cm radio luminosity to bolometric luminosity correlation, L ra( j — L\, a \, for the sample and 
discuss the nature of the radio emission in terms of the available sources of ionized gas. We 
also investigate the L raa - — Lir correlation and suggest that radio flux density may be used as a 
proxy for the internal luminosity of low luminosity protostars. 

Key words: Radiation mechanisms:general - ISM:general - ISMxlouds - stars:formation 



1 INTRODUCTION lengths of 1-2 cm (Drain & Lazarian 1998) and this is the wave- 

length range where most spinning dust observations are made. Al- 



The radio spectrum of dark clouds has recently come under scrutiny 
following the proposal that their centimetre-wave emission is due to 
spinning dust (Casassus et al. 2006; AMI Consortium: Scaife et al. 
2009a, b; 2010). These observational studies have largely concen- 



though a few observations towards specific sources have been made 
at 2cm (Rodriguez & Canto 1983; Bieging, Cohen & Schwartz 
1984; Pravdo 1985; Rodriguez & Reipurth 1997; Andre, Motte & 

Bacmann 1999) their number is still limited, particularly in the case 
trated on the arcminute scale emission, neglecting the known small „, , , . 

, .„ _ , of low luminosity protostars and Class O objects, 

scale radio emission detected from a significant fraction of pro- 
tostars (Anglada 1995; Rodriguez et al. 1989). As these investiga- Longer wavelength radio emission is capable of escaping from 
tions become more detailed it is necessary to quantify correctly this the high column density envelopes which surround deeply em- 
emission in order to avoid confusing it with that arising from spin- bedded young protostars, making it an appropriate mechanism for 
ning dust when observing at resolutions where the two may not be searching for such objects. Very high sensitivity observations are 
easily distinguishable. Most radio protostar searches and surveys required as the thermal (vibrational) dust spectrum falls off steeply 
(e.g. Anglada 1995; Stamatellos et al. 2009; Andre et al. 1987) are m intensity at longer wavelengths. However, young and low lumi- 
conducted at 3.6 and 6 cm wavelengths, however the spectrum of nos i ty protostars are also known to produce additional radio emis- 
emission due to spinning dust is thought to peak at shorter wave- s i on which can make them more easily detectable. In general this 

radio emission has been observed to possess a rising spectrum, in- 
dicating that it occurs as a consequence of free-free radiation from 

* We request that any reference to this paper cites "AMI Consortium: Scaife ionized § as - The s P ectral index of such emission (defined as a, 

etal. 2010" where S v °= V a ) can lie anywhere in the range —0.1 ^ a ^ 2. 

t E-mail: ascaife@cp.dias.ie An index of a = —0.1 indicates optically thin free-free emis- 
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sion, whereas a = 2 indicates optically thick; values intermediate 
to these limits imply a medium which is partially optically thick 
(Reynolds 1986). In the case of a spherically symmetric shell of 
decreasing density with n e <* r~ 2 , a ~ 0.6. A number of mech- 
anisms have been proposed to produce free-free emission in the 
immediate vicinity of stellar objects: where a high enough ionizing 
flux is present, generally in later type T Tauri stars, photoionization 
can support an embedded Hll region (Churchwell 1990); a fully 
ionized stellar wind, again associated more often with later type 
stars, would produce an easily observable radio signal (Panagia & 
Felli 1975; Curiel et al. 1989), as could a partially ionized, colli- 
mated outflow (Reynolds 1986); the accretion shock on the surface 
of protostellar discs may heat and ionize infalling gas (Winkler & 
Newman 1980; Cassen & Moosman 1981); the molecular outflows 
from young, low mass protostars and the neutral winds thought to 
initiate them can generate free-free emission through shock ion- 
ization as they impact on the surrounding envelope (Curiel et al. 
1987; 1989; Rodriguez & Reipurth 1996). In the case of low lumi- 
nosity protostars it is the last of these mechanisms which is consid- 
ered to be most plausible and indeed, where measured, the outflow 
force of most protostellar jets has been found to be energetically 
viable to explain the observed cm-wave radio emission (Anglada 
1995). This theory is also supported by very high angular reso- 
lution observations of the radio emission from protostars (Simon 
et al. 1983; Anglada 1996), which find that the radio emission is 
elongated along the direction of the molecular outflows. 

Recent measurements of the luminosity distribution of proto- 
stars from the Spitzer Space Telescope (Dunham et al. 2008; Evans 
et al. 2009) have aggravated the "luminosity problem" first artic- 
ulated by Kenyon et al. (1990). The trend is increasing numbers 
of protostars at lower luminosities down to an internal luminosity, 
which describes the total luminosity of the central protostar and 
circumstellar disc, of L; nt ~ 0.1 L©, below which the numbers start 
to decline. For young embedded protostars the accretion luminos- 
ity will dominate over that from the stellar photosphere. However, 
simple star formation models, such as the spherical accretion model 
of Shu (1977), predict that a low mass source on the stellar/brown 
dwarf boundary (M ~ 0.08 M©) should have an internal luminos- 
ity Li M — 1-6L© from accretion alone (Evans et al. 2009). Non- 
steady accretion, starting in the earliest protostellar stages, is cur- 
rently the best solution to this discrepancy (Kenyon & Hartmann 
1995; Young & Evans 2005; Enoch et al. 2007). The luminosity 
problem is most difficult to rectify in very low luminosity objects 
(VeLLOs; Young et al. 2004; Dunham et al. 2008) with extreme 
luminosities L mt ^ 0.1L©. The nature of these objects is unclear, 
whether they are young Class O protostars which are just powering 
up, or are more evolved but in a low accretion state (Dunham et al. 
2008; Evans et al. 2009). 

The VeLLO IRAM 04 191 -IRS (Andre et al. 1999) is a good 
example of this situation. IRAM 04191-IRS has a well measured 
molecular outflow, from which the mass accretion onto the proto- 
star can be calculated (Andre et al. 1999; Dunham et al. 2008). The 
luminosity expected from this accretion, using steady accretion ar- 
guments, is found to be 25 times the measured value (Dunham et al. 
2006). However, it is not the case that all VeLLOs are so well de- 
fined, and indeed many of their physical properties are often incon- 
sistent (Bourke et al 2006). The existing sample of VeLLO sources, 
although expanded by Spitzer (Dunham et al. 2008), is in no way 
complete. VeLLOs are difficult to confirm as protostars in the infra- 
red due to their low luminosity and embedded nature, and measur- 
ing their molecular outflows may also be problematic as exempli- 
fied by the case of L1014 (Bourke et al. 2005; Crapsi et al. 2005). It 



is also the case that VeLLOs are often found in cores which are not 
only assumed to be starless, but which were also not believed to be 
approaching collapse (Bourke et al. 2006). Nevertheless, identify- 
ing a complete sample of these low luminosty embedded protostars 
is vital for understanding low mass star formation. 

This paper presents high resolution, high sensitivity observa- 
tions towards a sample of known or proposed protostellar objects 
all located within clouds targeted by a joint Arcminute Microkelvin 
Imager (AMI) and Sunyaev-Zel'dovich Array (SZA; Carlstrom 
et al. ***) spinning dust dark cloud sample (Scaife et al. in prep) 
which was selected from the Spitzer small clouds and cores pro- 
gramme (Evans et al. 2003) and is designed to measure the emis- 
sion from the clouds on arcminute scales. The majority of the cores 
observed here have low bolometric luminosity, with a significant 
proportion being Very Low Luminosity Objects. Since these studies 
are on differing scales but target the same regions we now explic- 
itly define our nomenclature. Cores are defined as compact objects 
on scales of a few tenths of a parsec and are typically identified at 
sub-mm wavelengths. Cores can be protostellar or starless, a dis- 
tinction made on the basis of their SEDs, the presence of an IR 
source, a molecular outflow or a compact cm-wave source (Andre 
et al. 2000). VeLLOs are a subset of cores, which until recently 
were often presumed to be starless but are in fact a low luminos- 
ity embedded population. Cores identified in the sub-mm are suf- 
fixed "SMM", whilst those identified in the infra-red are suffixed 
"IRS". Clouds or dark clouds are regions of high visual extinction 
on scales of a few parsecs. They contain large quantitities of molec- 
ular gas, which often harbours cores and active star formation. On 
larger scales of a few parsecs to a few tens of parsecs we refer 
to extended regions of presumably interconnected clouds as cloud 
complexes. 

The organization of this paper is as follows. In § 2 we describe 
the sample of targets to be observed, and in § 3 we describe the 
AMI Large Array (AMI-LA) telescope, the observations and the 
data reduction process. In § 4 we comment upon the results of the 
observations and compare them to predictions. Detailed notes on 
individual fields are given in § 5. The nature of the radio emission 
is discussed in § 6 and correlations with the IR properties of the 
sample are derived in § 6.3. Finally we present our conclusions in 
§7. 



2 SAMPLE 

The thirteen fields observed by the AMI-LA are listed in Table 1. 
These fields were selected from the full AMI-SZA spinning dust 
sample (Scaife et al. in prep) on the basis of being known to con- 
tain cores, or alternatively to contain possible embedded objects 
from the Spitzer catalog of Dunham et al. (2008). Where neces- 
sary we shall identify candidates from this catalog by their catalog 
number, i.e. [DCE08]-nnn. These candidates were ranked by Dun- 
ham et al. (2008) as belonging to one of six "groups", with Group 
1 being those most likely to be true embedded objects, and Group 
6 those least likely. Table 1 lists the fields in order of increasing 
Right Ascension with the associated cloud named in Column [1] . 
In Columns [2] and [3] it lists the pointing center observed for each 
field, and in Column [4] the association of the cloud to a cloud 
complex. In Columns [5-7] the existing sub-mm observations are 
listed, the references for which are given in Column [8]. Column 
[9] lists the number of cores of varying class thought to be in each 
field. Within the AMI-LA fields there are 21 known cores or po- 
tential embedded objects. Five fields contain known VeLLO ob- 
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jects: IRAM 04191 (Dunham et al. 2006), L1521F (Bourke et al. 
2006; Terebey et al. 2009; Shinnaga et al. 2009), L673-7 (Dun- 
ham et al. 2008), LI 148 (proposed; Kauffmann et al. 2005), and 
L1014 (Young et al. 2004; Bourke et al. 2005; Shirley et al. 2007). 
To this list we add two further potential VeLLO objects: L723- 
IRS ([DCE08]-026) and L1165-IRS ([DCE08]-039), which were 
included in the catalog of Dunham et al. (2008) but ranked as Group 
6. 

All of the fields have been at least partially observed in the 
sub-mm, although three of the VeLLO objects (L723-IRS and 
LI 165-IRS) are not covered by existing observations. Several fields 
also have data at lower radio frequencies: IRAM 04191 (Andre, 
Motte & Bacmann 1999), L1521F (Harvey et al. 2002), B35A 
(Rodriguez et al. 1989), L723 (Anglada et al. 1991; Carrasco- 
Gonzalez et al. 2008), L1014 (Shirley et al. 2007) and L1221 
(Rodriguez & Reipurth 1998; Young et al. 2009). The correspon- 
dence of the AMI-LA measurements to these observations are dis- 
cussed in more detail in § 6.2. 



3 OBSERVATIONS 

AMI comprises two synthesis arrays, one of ten 3.7 m antennas 
(Small Array; SA) and one of eight 13 m antennas (Large Ar- 
ray; LA), both sited at Lord's Bridge, Cambridge (AMI Consor- 
tium: Zwart et al. 2008). The telescope observes in the band 13.5- 
17.9 GHz with cryostatically cooled NRAO indium-phosphide 
front-end amplifiers. The overall system temperature is approxi- 
mately 25 K. Amplification, equalization, path compensation and 
automatic gain control are then applied to the IF signal. The back- 
end has an analogue lag correlator with 16 independent correla- 
tions formed at twice Nyquist rate for each baseline using path de- 
lays spaced by 25 mm. From these, eight complex visibilities are 
formed in 0.75 GHz bandwidth channels. In practice, the two low- 
est frequency channels (1 & 2) are not generally used due to a lower 
response in this frequency range and interference from geostation- 
ary satellites. The data in this paper were taken with the AMI Large 
Array. 

Observations of the thirteen fields listed in Table 1 were made 
with the AMI-LA between October 2009 and July 2010. The clouds 
were observed as single pointings, with the exception of B35A and 
LI 148. B35A was observed as a 7 field mosaic in order to cover 
the same sky area as earlier VLA observations, see § 5, and LI 148 
was observed in two separate pointings as the two embedded cores 
proposed by Dunham et al. (2008) for this cloud are separated by 
more than the AMI-LA primary beam FWHM. 

AMI-LA data reduction is performed using the local software 
tool REDUCE. This applies both automatic and manual flags for 
interference, shadowing and hardware errors, Fourier transforms 
the correlator data to synthesize frequency channels and performs 
phase and amplitude calibrations before output to disc in uv FITS 
format suitable for imaging in MPS. Flux (primary) calibration is 
performed using short observations of 3C286 and 3C48. We as- 
sume I+Q flux densities for this source in the AMI LA channels 
consistent with the updated VLA calibration scale (Rick Perley, pri- 
vate comm.), see Table 2. Since the AMI-LA measures I+Q, these 
flux densities include corrections for the polarization of the calibra- 
tor sources. A correction is also made for the changing interven- 
ing air mass over the observation. From other measurements, we 
find the flux calibration is accurate to better than 5 per cent (AMI 
Consortium: Scaife et al. 2008; AMI Consortium: Hurley-Walker 
et al. 2009). Additional phase (secondary) calibration is done us- 



Table 2. AMI-LA frequency channels and primary calibrator flux densities 
measured in Jy. 



Channel No. 


3 


4 


5 


6 


7 


8 


Freq. [GHz] 


14.27 


14.99 


15.71 


16.43 


17.15 


17.87 


3C48 


1.85 


1.75 


1.66 


1.58 


1.50 


1.43 


3C286 


3.60 


3.54 


3.42 


3.31 


3.21 


3.11 


3C147 


2.75 


2.62 


2.50 


2.40 


2.30 


2.20 



ing interleaved observations of calibrators selected from the Jodrell 
Bank VLA Survey (JVAS; Patnaik et al. 1992). After calibration, 
the phase is generally stable to 5° for channels 4-7, and 10° for 
channels 3 and 8. The FWHM of the primary beam of the AMI LA 
is ~ 6 arcmin at 16 GHz. 

Reduced data were imaged using the AIPS data package. 
CLEAN deconvolution was performed using the task IMAGR which 
applies a differential primary beam correction to the individual fre- 
quency channels to produce the combined frequency image. The 
broad spectral coverage of AMI allows a representation of the spec- 
trum between 14.3 and 17.9 GHz to be made independently of other 
telescopes when sufficient signal to noise is present and in what 
follows we use the convention: S v °= v a , where S v is flux density 
(rather than flux, F v = v5 v ), V is frequency and a is the spectral 
index. All errors quoted are 1 o~. 

The observations towards the thirteen fields listed in Table 1 
and described in § 2 are summarized in Table 3. The details of 
each individual observation including the date and both the pri- 
mary and secondary calibration sources are listed, along with the 
resulting rms noise level in the map and the dimensions of the nat- 
urally weighted synthesized beam. The rms noise level varies be- 
tween fields due to the different levels of data flagging required 
following periods of poor weather conditions or interference from 
non-astronomical sources, such as geostationary satellites. Where 
possible for observations that were heavily flagged, a second obser- 
vation was made and these are also indicated in Table 3. 



4 RESULTS 

Within the 13 fields observed by the AMI-LA we detect 40 sources, 
and we identify 15 of these as being associated with the 21 possible 
cores in these areas. The combined channel maps from the AMI- 
LA are shown in Figs. 1-5. The rms noise on each map, (Tj-jj^ , is 
found using the AIPS task IMEAN which fits a Gaussian profile to 
the histogrammed pixel values for each field. These values may be 
found in Table 3 and in the figure captions. We identify all objects 
with peak flux densities greater than 5o" rms , where that peak lies 
within the FWHM of the AMI-LA primary beam, as being true 
sources and these are listed in Table 4. The two detection conditions 
are relaxed for some individual cases and where this occurs it is 
indicated in Table 4 and commented upon in Section 5. 

Table 4 lists sources in order of increasing Right Ascension for 
each field, with peak positions given in Columns [3] and [4]. Given 
the low signal-to-noise in the individual channel data for these tar- 
gets, the reduced phase stability of AMI-LA Channel 3 was found 
to produce an rms noise too high relative to the other channels to 
contribute usefully to the combined dataset, consequently it has not 
been used in this work. In addition, since channel 8 of the AMI-LA 
is prone to satellite interference it is sometimes not used in mak- 
ing the final combined bandwidth images. As a rule, this channel is 
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Table 1. AMI-LA Fields. Column [1] Name of field, [2] Right Ascension, [3] Declination, [4] Association, [5] flag on existing data from SCUBA, [6] flag on 
existing data from MAMBO, [7] flag on existing data from SHARC-II, [8] references for [5-7], [9] number of cores (confirmed and proposed) within the area 
covered by the AMI-LA observation. 



Name 


RA 


Dec 


Association 


SCUBA 


MAMBO 


SHARC-II 


refs. 


Ncore 




(12000) 


(J2000) 




Y/N 


Y/N 


Y/N 






IRAM04191 


04 21 57 


15 29 46 


Taurus 


Y 


N 


Y 


1,6 


2 


L1521F 


04 28 39 


26 51 36 


Taurus 


Y 


Y 


Y 


1,5,6 


1 


B35A 


05 44 29 


09 08 57 




Y 


N 


Y 


1,5,6 


2 


L673 


19 20 26 


11 22 19 


Aquila Rift, Cloud B 


Y 


N 


N 


2 


3 


CB188 


19 20 15 


11 36 08 


Aquila Rift, Cloud B 


N 


Y 


Y 


4,5,6 


1 


L673-7 


19 21 35 


11 21 23 


Aquila Rift, Cloud B 


N 


Y 


Y 


5,6 


1 


L723 


19 17 44 


19 15 24 




Y 


N 


Y 


1,6 


2 


L1152 


20 35 46 


67 53 02 




Y 


N 


Y 


1,6 


1 


L1148 


20 40 57 


67 23 05 


Cepheus Flare 


N 


Y 


Y 


3,5,6 


2 


BERN 48 


20 59 15 


78 22 60 


Cepheus Flare 


Y 


Y 


Y 


5,6 


1 


L1014 


21 24 08 


49 59 09 




Y 


Y 


Y 


2,5,6 


1 


LI 165 


22 06 50 


59 02 46 


Cloud 157, Cyg OB7 


N 


N 


Y 


2,6 


2 


L1221 


22 28 07 


69 00 39 




Y 


N 


Y 


1,6 


2 



[1] Young et al. 2006; [2] Visser et al. 2002; [3] Kirk, Ward-Thompson & Andre 2007; [4] Launhardt et al. 1997; [5] Young et al. 2006b; [6] Wu et al. 2007. 



Table 3. AMI-LA Observations. Column [1] Name of field, [2] Date of initial observation, [3] Primary Calibrator, [4] Secondary Calibrator", [5] Date of 
additional observation, [6] Primary Calibrator, [7] Secondary Calibrator , [8] AMI-LA combined data synthesized beam FWHM major axis, [9] AMI-LA 
combined data synthesized beam FWHM minor axis, and [10] rms noise fluctuations on the combined channel map. 



Name 


Date 
(dd-mm-yy) 


1° 


2° 


Date 
(dd-mm-yy) 


1° 


2° 


A0 ma j 
(arcsec) 


A0 ra in 

(arcsec) 


^rms 
bm ' 


IRAM04191 


09-10-09 


3C48 


J0424+1442 


07-11-09 


3C286 


J0424+1442 


49.9 


27.4 


17 


L1521F 


17-10-09 


3C147 


J0426+2350 








42.0 


26.3 


16 


B35A* 


08-11-09 


3C48 


J055 1+0829 








53.6 


28.0 


67 


L673 


13-12-09 


3C48 


J1922+1530 








36.5 


31.2 


31 


CB188 


27-12-09 


3C48 


J1922+1530 








46.0 


21.9 


24 


L673-7 


18-10-09 


3C286 


J1922+1530 








32.0 


26.3 


16 


L723 


14-12-09 


3C286 


J1905+1943 


08-10-09 


3C286 


J1925+2106 


42.9 


24.0 


29 


L1152 


26-12-09 


3C286 


J2035+5821 


05-01-10 


3C286 


J2035+5821 


38.8 


22.6 


49 


L1148 c 


12-10-09 


3C286 


J2052+6858 


20-10-09 


3C286 


J2052+6858 


36.6 


26.5 


22 


BERN 48 


15-11-09 


3C48 


J2051+7441 








30.0 


27.4 


23 


L1014 


14-10-09 


3C48 


J2123+4614 








31.2 


25.4 


21 


LI 165 


30-06-10 


3C147 


J2223+6249 


02-07-10 


3C48 


J2223+6249 


32.3 


29.0 


21 


L1221 


11-10-09 


3C286 


J2230+6946 








27.4 


25.1 


19 



" Secondary calibrators are selected from the JVAS catalog (Patnaik et al. 1992). 

* The 7 pointings of the B35A cloud mosaic were observed in an interleaved fashion during single run. 

c The 2 pointings towards the LI 148 cloud were observed on separate dates, the later observation corresponds to the southern pointing. 



only included when the resulting noise level in the combined chan- 
nel map, O^ 8 , is lower than that of a map using only channels 4 
to 7, . As this varies from field to field, changing the nomi- 
nal observation frequency slightly, Column [5] lists this frequency 
for each dataset. For ease of notation AMI-LA measurements will 
henceforth be denoted as being made at 16 GHz (1.8 cm) when de- 
scribed in the text; for the more precise measurement frequency of 
individual fields we refer the reader to Table 4. Column [6] lists 
the primary beam corrected peak flux density for each detected 
source with no correction for the local background, and Column 
[7] the primary beam corrected integrated flux density. Integrated 
flux densities are found using the FITFLUX program (Green 2007). 
This method calculates flux densities by removing a tilted plane 
fitted to the local background and integrating the remaining flux 
density. This is done by drawing a polygon around the source and 
fitting a tilted plane to the pixels around the edge of the polygon. 
Where an edge of the polygon crosses a region confused by an- 



other source, the background is subjective and this edge is omitted 
from the fitting. Since the extracted flux density is dependent to 
some degree on the background emission, we repeat this process 
using five irregular polygons, each varying slightly in shape. The 
final flux density is the average of that extracted from these five 
polygons. The error on the integrated flux densities is calculated as 

(7s = \j (0.055; nt ) 2 + tf| t + (7^ ms , where <7fi t is the standard devia- 
tion of the flux densities found from the five polygonal apertures, 
(7 rms is the rms noise determined using IMEAN and 0.05Sj nt is a con- 
servative 5 per cent absolute calibration error. Where the 5 per cent 
calibration uncertainty is not dominant, errors determined for this 
sample are in general heavily dependent on the rms noise, C7 rms , 
rather than the fitting error, <7 fit , which is found to be small. Col- 
umn [8] of Table 4 then identifies which of the detected sources are 
associated with known cores. 
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4.1 A Note on Expected Sources 

At 16 GHz we expect a certain number of extragalactic radio 
sources to be seen within each of our fields. To quantify this 
number we use the 15 GHz source counts model from de Zotti 
et al. (2005) scaled to the extended 9C survey source counts (Wal- 
dram et al. 2010). The average rms noise from our datasets is 
~ 25/iJy beam -1 and from this model we predict that we should 
see 0.07 sources arcmin , ~ 2 radio sources within a 6arcmin 
FWHM primary beam above a 5 o~ flux density of 1 25 fl Jy, or 3 1 ± 6 
sources in total. On average in our sample we see 0.09 sources 
arcmin -2 . Since we identify 15 of our 40 sources with protostellar 
cores, not accounted for by the extragalactic radio source counts 
model, this is consistent with the discrepancy in the measured and 
predicted source counts. Since the short frequency coverage of the 
AMI-LA band makes spectral index estimates highly uncertain for 
low signal to noise detections, and given this statistic, we tenta- 
tively identify any sources found in the target fields that do not 
have a known protostellar association as extragalactic. 

5 NOTES ON INDIVIDUAL FIELDS 
5.1 IRAM04191 

IRAM 04191+1522 (IRAM 04191-IRS), see Fig. 1 (a), was iden- 
tified as a Class O source by Andre, Motte & Bacmann (1999; 
AMB99). It has a large extended molecular outflow, seen in CO(2- 
1) (AMB99) and a flattened core mapped in N2H+ (Dunham et al. 
2006). AMB99 reported radio flux densities for this object at 5, 
8 and 15 GHz from the VLA, which they stated were consistent 
with a spectral index of a = 0.6, similar to that seen in the radio 
from VLA1623. We find a flux density of 5 i6 = 0. 173 ± 0.020/iJy 
for IRAM 04191-IRS, consistent with the earlier VLA measure- 
ment of S| 5 LA = 0. 16 ± 0.05 AtJy- Combining the VLA and AMI- 
LA channel data gives a spectral index of a\\ 5 = 0.45 ± 0.20. 
Across the AMI-LA band alone we measure a spectral index of 
alJj = 1.44 ± 1.09. The large error from the relatively small fre- 
quency coverage means that this is broadly consistent with that 
found at lower frequencies, but it may also indicate a steepening 
of the spectrum. The morphology of the object in the AMI-LA 
maps has a deconvolved major axis of 29" with a position angle 
of 1 12° and is unresolved in the minor direction. This morphology 
is similar to that seen at 350 flm (Dunham et al. 2006). Also within 
the AMI-LA field, to the north-east of IRAM 04191-IRS, is the 
similarly named IRAS 04191, a Class I source. We detect no radio 
emission from this object. There is only one further source detected 
above 5 a within the AMI-LA primary beam FWHM towards this 
field, and in the absence of further information we identify this as 
an extragalactic radio source. 



5.2 L1521F 

L1521F-IRS, see Fig. 1 (b), is thought to be a core in the very 
earliest stages of gravitational collapse (Terebey et al. 2009). Un- 
like the VeLLO sources IRAM04191-IRS or L1014-IRS it does 
not possess a detected molecular outflow. Previous deep cm-wave 
observations of the L1521F cloud (Harvey et al. 2002), which 
covered L1521F-IRS, at 3.6 cm wavelength did not find a coin- 
cident radio source and instead used a nearby object, which we 
identify as NVSS 042841+265355, to place an upper limit on 
the bolometric luminosity of the core. In Fig. 1 (b) the source 



NVSS 042841+265355 can be seen to the north of the phase cen- 
tre, just within the primary beam FWHM. This source has a primary 
beam corrected flux density of = 6. 1 ± 0.3 mjy and in combi- 
nation with Si.4 = 37.3 ± 1.2mJy (NVSS) and 5 8 . 35 = 10.8mJy 
(Harvey et al. 2002) we find this consistent with the presence of 
an extra-galactic non-thermal point source (as also concluded by 
Harvey et al). However, in addition to NVSS 042841+265355 we 
also detect a two component radio structure more coincident with 
L1521F-IRS. Shinnaga et al. (2009) detected extended warm gas 
(30-70 K) towards L1521F-IRS and have suggested that it is highly 
likely this gas is heated through shocks produced by an outflow- 
ing/rotating component impacting on the surrounding cold dynam- 
ically collapsing medium. Such shocks would also be likely to pro- 
duce radio emission from shock ionization and it may be that the 
double radio structure we see here is a consequence of L1521F- 
IRS's bipolar outflow, although we note that the radio emission is 
not aligned east-west following the morphology of the outflow as 
mapped by Bourke et al. (2006). This may be a consequence of 
the non-uniformity of the surrounding interstellar medium, or as 
suggested by Shinnaga et al. (2009) be due to the outflows hav- 
ing rotated from the position of the shocks responsible for the radio 
emission. Of the three further sources in this field, none have known 
radio or protostellar counterparts and in the absence of further in- 
formation we identify them as extragalactic sources. 



5.3 B35A 

B35A SMM-1, see Fig. 2, is an evolved protostar with observed 
molecular outflows (Myers et al. 1988). A 6 cm radio survey with 
the VLA (Rodriguez et al. 1998) placed an upper limit of 0.2 mjy 
flux density on any radio emission from the protostar and we 
here observe a 7-field mosaic of the same sky area. At 16 GHz 
we observe 0.94 ± 0.08 mjy of flux density coincident with the 
B35A SMM-1 core, as identified by the Submillimeter High An- 
gular Resolution Camera II (SHARC-II; Wu et al. 2007). Al- 
though the radio source extends slightly in their direction we do 
not see separate peaks from the adjacent cores B35A SMM-2 and 
B35A SMM-3. We detect all the radio sources identified at 6 cm 
with the VLA and confirm that they are universally steep-spectrum 
objects and therefore most likely to be non-thermal extragalactic 
sources. The second of the candidate protostars proposed by Dun- 
ham et al. (2008), which we denote B35A-IRS, is located in the 
most south-easterly AMI-LA pointing adjacent to a bright extra- 
galactic radio source. Any emission from this source cannot be 
separately resolved by the AMI-LA synthesized beam. Of the three 
further sources detected in ths field, none have known associations 
and we identify them as extragalactic sources. 



5.4 L673 

L673, see Fig. 1 (c), is a large star forming region in the Aquila Rift 
and the field observed here covers only part of it. Radio counter- 
parts are detected for the known sub-mm cores L673 SMM-1 and 
L673 SMM-2 (Visser et al. 2002). A third sub-mm core, identified 
as L673 SMM-7 in Visser et al. (2002), which lies to the north-west 
of L673 SMM-1 is not detected. This core has no detected molecu- 
lar outflow (Visser et al. 2002) and is assumed starless. The two ra- 
dio counterparts detected in the AMI-LA maps have been denoted 
L673 SMM-1 and L673 SMM-2 in Table 4. However we note that 
although L673 SMM-1 is coincident with the sub-mm core detected 
in the SCUBA data of Visser et al. (2002), the radio counterpart to 
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Table 4. Sources detected within the AMI-LA fields. 



Field 


Source 


RA 
(J2000) 


Dec 
(J2000) 


V 

(GHz) 


S eak 
(jtiJybeairr 1 ) 


(j"Jy) 


assoc. 


IRAM04191 


(1) 


04 21 56.7 


15 29 36.0 


16.43 


164.3 


172.9±20.1 


IRAM 04191-IRS 




(2) 


04 22 05.0 


15 28 01.0 


16.43 


185.5 


316.7±27.1 




L1521F 


(1) 


04 28 36.4 


26 53 06.0 


16.07 


120.0 


90.5±17.9 






(2) 


04 28 38.6 


26 51 11.0 


16.07 


132.9 


139.2±18.2 


L1521F-IRS 




(3) 


04 28 40.9 


26 51 51.0 


16.07 


138.8 


234.2±22.0 


L1521F-IRS 




(4) 


04 28 41.2 


26 53 56.0 


16.07 


5822.6 


6136.1±308.4 






(5) 


04 28 45.4 


26 53 56.0 


16.07 


247.4 


289.8±19.9 






(6) 


04 28 47.2 


26 52 36.0 


16.07 


117.5 


139.1±18.7 




B35A 


(1*) 


05 44 06.0 


09 08 42.0 


16.07 


1110.8 


936.1 ±82.1 






(2*) 


05 44 08.3 


09 08 02.0 


16.07 


1440.3 


1433.5±100.1 






(3) 


05 44 19.5 


09 11 17.0 


16.07 


352.3 


363.7±70.2 






(4) 


05 44 22.5 


09 09 47.0 


16.07 


585.6 


558.2±73.1 






(5) 


05 44 29.6 


09 08 52.0 


16.07 


591.7 


943.0±84.5 


B35A SMM-1 




(6) 


05 44 31.0 


09 07 42.0 


16.07 


355.7 


683.2±78.4 






(7) 


05 44 43.8 


09 03 52.0 


16.07 


811.3 


1139.2±92.3 






(8) 


05 44 48.5 


09 03 27.0 


16.07 


7450.2 


8007.4±4()7.3 






(9) 


05 44 53.2 


09 06 22.0 


16.07 


676.7 


642.3±75.4 






(10) 


05 44 53.9 


09 10 17.0 


16.07 


389.2 


430.9±76.4 






(11) 


05 45 02.0 


09 05 01.9 


16.07 


394.1 


343.7±70.1 




L673 


(1) 


19 20 25.0 


11 22 15.0 


16.07 


354.6 


455.2±39.3 


L673 SMM-1 




(2) 


19 20 26.0 


11 20 15.0 


16.07 


218.4 


289.7±35.9 


L673 SMM-2 




(3) 


19 20 35.2 


11 21 00.0 


16.07 


300.6 


380.3±39.0 




CB188 


(1*) 


19 20 14.7 


11 35 48.0 


16.43 


179.5 


179.5* 


CB188 SMM-1 


L673-7 


(1) 


19 21 28.9 


11 22 48.0 


16.07 


163.6 


548.1±37.5 






(2) 


19 21 33.0 


11 22 48.0 


16.07 


184.6 


428.4±34.3 




L723 


(1) 


19 17 44.4 


19 15 24.0 


16.43 


122.2 


220.0±32.() 


L723-IRS 




(2) 


19 17 46.5 


19 14 29.0 


16.43 


127.3 


154.0±30.8 






(3*) 


19 17 53.5 


19 12 29.0 


16.43 


610.0 


689.5±45.5 


L723 SMM-1 


L1152 


(1) 


20 35 29.2 


67 52 07.0 


16.43 


2470.2 


2471.5±133.1 




LI 148 


(1) 


20 40 41.4 


67 25 20.0 


16.07 


116.7 


205.5±26.3 






(2+) 


20 40 57.0 


67 22 55.0 


16.07 


83.7 


71.9+ 


L1148-IRS 


BERN 48 


(1) 


20 59 15.0 


78 23 04.9 


16.43 


597.4 


615.2±38.7 


Bern 48 


L1014 


(1) 


21 23 55.0 


49 58 59.0 


16.07 


1615.6 


1680.2±87.1 






(2) 


21 23 56.1 


49 57 54.0 


16.07 


203.3 


249.7±28.7 






(3) 


21 24 08.5 


49 59 04.0 


16.07 


166.3 


298.5±26.5 


L1014-IRS 




(4) 


21 24 17.3 


50 01 04.0 


16.07 


242.7 


229.8±28.8 






(5) 


21 24 20.4 


49 57 24.0 


16.07 


682.3 


634.0±38.7 




LI 165 


(1) 


22 06 17.5 


59 04 29.9 


16.43 


250.5 


397.1 ±33.0 






(2) 


22 06 30.5 


59 02 40.0 


16.43 


109.9 


81.6±22.8 






(3) 


22 06 35.0 


59 02 35.0 


16.43 


186.4 


346.4±30.6 


L1165-IRS 




(4) 


22 06 49.9 


59 03 00.0 


16.43 


240.7 


305.8±28.2 


LI 165 SMM-1 


L1221 


(1) 


22 27 43.7 


69 00 48.9 


16.43 


1535.5 


1704.3±87.7 






(2) 


22 27 51.2 


69 00 24.0 


16.43 


187.9 


182.2±24.2 






(3) 


22 27 52.1 


69 01 39.0 


16.43 


1867.5 


2081.7±106.4 






(4) 


22 28 03.3 


69 01 14.0 


16.43 


282.8 


358.8±27.6 


L1221-IRS1 




(5) 


22 28 07.0 


69 00 39.0 


16.43 


367.4 


358.9±28.1 


L1221-IRS3 




(6) 


22 28 35.8 


69 00 08.9 


16.43 


350.9 


353.9±27.8 





*These sources are outside the AMI-LA primary beam FWHM. 

tThese sources are detected at 3 o". 

* Lower limit on flux density, see text for details. 



L673 SMM-2 is offset slightly to the north of the sub-mm posi- 
tion. The peak flux density of the L673 SMM-1 radio counterpart 
is also coincident with the candidate embedded object [DCE08]- 
027, which is a Group 1 candidate according to the classification 
of Dunham et al. (2008). A second candidate [DCE08]-028 (Group 
5) is offset very slightly to the north-east, away from the peak. Al- 
though the sub-mm core L673 SMM-2 is very close in position to 
the candidate embedded object [DCE08]-029 (Group 1) the slight 
offset of the radio emission to the north makes the radio counter- 
part more coincident with [DCE08]-030 (Group 3). However the 



two candidate sources are separated by less than one synthesized 
beam and it is difficult to correctly assign the radio emission to a 
single candidate with observations of this resolution. 



Aside from the previously mentioned sources there is one fur- 
ther object in this field, which has no known radio or protostel- 
lar counterpart and therefore we ascribe to being an extragalactic 
source. 
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Figure 1. (a) IRAM 04191; (b) L1521F; (c) L673; and (d) CB188. AMI-LA data are shown as greyscale and contours linearly from 3 cr, where a = 17, 16,26 
and 24/jJy beam -1 respectively. Spitzer positions and candidate embedded objects from Dunham et al. (2008) are indicated by crosses, and the positions of 
sub-mm cores are shown as unfilled squares, see text for details. The AMI-LA primary beam FWHM is shown as a circle and the AMI-LA synthesized beam 
is shown as a filled ellipse in the lower left corner of each map. 



5.5 CB188 

CB188, see Fig. 1 (d), shows an elongated radio structure with a 
slight localized enhancement at the position of the sub-mm core 
(CB 1 88 SMM) as detected by the Max Planck Millimetre Bolome- 
ter (MAMBO; Kauffmann et al. 2008) and SHARC-II (Wu et al. 
2007). CB188 SMM-1 is a known Class I source (Launhardt et al. 
2010) and has molecular outflows which were measured by Yun & 
Clemens (1994). The radio emission seen at 16 GHz by the AMI- 
LA shows the same spur of emission extending to the north away 
from the position of the sub-mm core, which is seen at 450, 850 
and 1300 flm (Launhardt et al. 2010). However the bulk of the ex- 
tended radio emission seen in the AMI-LA map of this object is ex- 
tended East-West, whereas the molecular outflows associated with 
CB188 SMM-1 are seen to be almost exactly coincident with the 
core and therefore are thought to be "pole-on", and to extend along 
the line of sight towards this object (Yun & Clemens 1994). It is 
consequently difficult to associate the localized enhancements of 



radio emission along this extension as being associated with the 
outflows. Given the unusual extended nature of this emission com- 
pared with the other fields in our sample we are hesitant to associate 
it entirely with CB 1 88 SMM. Taking this into account, the complex 
radio morphology of this object prevents us from extracting an in- 
tegrated flux density for CB188 SMM-1 alone, and so instead we 
place a lower limit on the integrated radio flux density, correspond- 
ing to the peak flux density at the position of the core. 

5.6 L673-7 

Although L673-7-IRS, see Fig. 3 (a), has previously been assumed 
to be a starless core, Dunham et al. (2008) classified it as a Group 
1 protostellar candidate based on its Spitzer data, albeit a low lumi- 
nosity Lir = 0.017 L© one. In spite of this we detect no source at 
the position of L673-7-IRS ([DCE08]-031). The two sources which 
are detected have no known associations and we tentatively identify 
them as extragalactic. 
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Figure 2. Primary beam corrected mosaic of B35A. AMI-LA data are shown as greyscale and contours linearly from 3 a, where a = 67 /iJy beam~ 1 , to 
lmJybeanT 1 and then incrementally in units of lmJybeanT 1 . Spitzer candidate embedded objects from Dunham et al. (2008) are indicated by vertical 
crosses (+), and the positions of sub-mm cores are shown as unfilled squares, see text for details. The location of B35A SMM-1 is highlighted by a square 
frame to the north-west of the centre. Radio sources identified by RR98 are marked by rotated crosses (x). The AMI-LA primary beam FWHMs are shown as 
circles and the AMI-LA synthesized beam is shown as a filled ellipse in the lower left corner. 



5.7 L723 

L723 SMM-1, see Fig. 3 (b), is a known Class O object with a very 
powerful molecular outflow (Bontemps et al. 1996). As its bolo- 
metric luminosity is also quite low, this places it in the early stages 
of Class O evolution. We detect L723 SMM-1 just outside the AMI- 
LA primary beam FWHM coincident with its known sub-mm posi- 
tion. The AMI-LA field is centred on [DCE08]-026, which we shall 
denote L723-IRS, and which is classified as a Group 6 embedded 
protostellar candidate by Dunham et al. (2008). Group 6 candidates 
are known to be not associated with regions of high volume den- 
sity and are therefore considered to be unlikely candidates, with the 
caveat that most of the dust continuum surveys used to determine 
regions of high volume density are limited to M > 0.1 — I.OMq. 
At 16 GHz we see a clear radio association to L723-IRS, with a 
second object offset to the south which may also be associated. 

5.8 L1152 

We detect no source at the position of L1152-IRS, see Fig. 3 (c). 
Only one source is visible within the AMI-LA primary beam which 
we identify with NVSS 203529+675207. This source has a primary 
beam corrected flux density of Sig = 2.47 ± 0. 13 mjy, which im- 
plies a spectral index of a\ ^ = 0. 15 ± 0.03. Due to its position off- 
set from LI 152-IRS we do not associate it with the core. However, 
we note that observations in the direction of LI 152-IRS with the 
AMI-LA are heavily affected by radio interference and the noise 
level of this observation is high compared to the rest of the sample. 

5.9 L1148 

The LI 148 cloud appears as two filaments of emission in the 
MAMBO 1.2 mm dust maps of this region (Kauffmann et al. 2008). 



We make a weak detection (3rj) of some extended radio emis- 
sion coincident with the northern most filament, see Fig. 4, also re- 
ferred to as L1148B in SHARC-II observations of this region. This 
northern filament is also coincident with a bright Spitzer source at 
24 and 70 /xm, LI 148-IRS, which has been proposed as a VeLLO 
(Kauffmann et al. 2005). Unlike other VeLLOs, such as L1014- 
IRS and IRAM 04191-IRS, no molecular outflow has yet been de- 
tected from LI 148-IRS and its status as a sub-stellar mass Class O 
source is still debatable. Dunham et al. (2008) place it in Group 3 
of their classification scheme, making it a probable protostar. If the 
enhancement of radio emission seen by AMI-LA is associated with 
ionization from a previously undetected molecular outflow or stel- 
lar wind from LI 148-IRS this would support the case for its classi- 
fication as a VeLLO. The second of the Spitzer protostar candidates 
for the LI 148 cloud, classified as Group 6, shows no corresponding 
radio emission in the AMI-LA data. With no known association, we 
identify the one further source detected in this field as extragalactic. 



5.10 Bern 48 

Bern 48 (also known as RN0129 and HH198) is a borderline Class 
I/II object (Connelley et al. 2008; Wu et al. 2008). The AMI-LA 
observations towards Bern 48, see Fig. 3 (d), show a slightly re- 
solved compact source (f2 src ~ 5") with a rising spectrum. The 
source is clearly extended in the sub-mm and on larger scales 
displays a cometary morphology in the 1.2 mm MAMBO maps 
(Kauffmann et al. 2008). The AMI-LA measures a flux density of 
5i6 = 615 ±39 jllJy, with a spectral index measured across the AMI 
band, although poorly constrained, of OfJJn = 1.25 ±0.96. A spec- 
tral index of ~ 1 is consistent with that seen towards hypercompact 
Hll regions, or indeed which may be expected from a stellar wind. 
If Bern 48 is a HCHlI region of consistent density this would imply 
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Figure 3. (a) L673-7; (b) L723; (c) L1152; and (d) BERN 48. AMI-LA data are shown as greyscale and contours linearly from 3 ff, where a = 19,22,49 
and 23 iijy beam - 1 . Spitzer candidate embedded objects from Dunham et al. (2008) are indicated by crosses, and the positions of sub-mm cores are shown as 
unfilled squares, see text for details. The AMI-LA primary beam FWHM is shown as a circle and the AMI-LA synthesized beam is shown as a filled ellipse in 
the lower left corner. 



an emission measure of > 10 9 pccirr 6 and a mass of ^ 100 M Q . 
Mass estimates for Bern 48, derived from MAMBO flux densities 
(Kauffmann et al. 2008) find M tot = 2.06 ± 0.02 M e inconsistent 
with this prediction. 

5.11 L1014 

The sub-mm core in L1014 was revealed not to be starless, as 
previously thought (Visser et al. 2002), by the Spitzer c2d pro- 
gramme and became the first VeLLO object (L1014-IRS; Lir ~ 
0.09 L ; Young et al. 2004). Deep VLA observations at 6 and 
3.6 cm revealed a radio counterpart with a rising spectral index 
of afH = 0.37 ±0.34. The AMI-LA measures a flux density of 
5i6 = 299±27,uJy towards L1014-IRS, constraining the spectral 
index to be Q^gg = 0.71 ±0.11. We observe the emission from 
L1014-IRS to be slightly extended, see Fig. 5 (a), suggesting once 
more that the radio emission arises not from the central engine itself 
but from the more extended envelope or associated outflow. 



There are also a number of bright radio sources in this field 
which can be found in the 1.4 GHz NVSS catalogue and that we 
identify as steep spectrum radio sources. The two sources detected 
in this field which do not have an NVSS counterpart we also choose 
to identify as extragalactic as they have no known protostellar as- 
sociations. 

5.12 L1165 

L1165-IRS ([DCE08]-039) is classified as a Group 6 source by 
Dunham et al. (2008), making it an unlikely candidate for an em- 
bedded object. In spite of this we observe an extended region of 
radio emission with a double-lobe structure, aligned east-west. Al- 
though the peak of neither lobe is exactly coincident with LI 165- 
IRS we cautiously associate the larger lobe, which is nearer to 
the FIR source with this object. A second potential Group 6 core 
in LI 165 ([DCE08]-040), with very low IR-luminosity, L IR = 
0.008 L , is listed by Dunham et al. (2008) to the north of LI 165- 
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Figure 5. (a) L1014; (b) LI 165; and (c) L1221. AMI-LA data are shown as greyscale and contours. In (a) & (b) contours are linear from 3 o~, where 
a = 21 jtJybeam -1 in both maps. In (c) contours are at 3, 6, 12, 24o~, where o" = 19/iJybeam - '. Spitzer candidate embedded objects from Dunham et al. 
(2008) are indicated by crosses, and the positions of sub-mm cores are shown as unfilled squares, see text for details. The AMI-LA primary beam FWHM is 
shown as a circle and the AMI-LA synthesized beam is shown as a filled ellipse in the lower left corner of each map. 



IRS and we see no emission towards this object. The two detected 
sources which we do not associate with L1165-IRS do have any 
known radio or protostellar counterparts and we therefore identify 
them as extragalactic sources. 

LI 165 SMM-1, to the east of L1165-IRS, was identified as a 
Class I object by Visser et al. (2002). The AMI-LA measurements 
of this source, see Fig. 5 (b), show a relatively compact, isolated 
source slightly extended in the north-south direction. The spectrum 
of LI 165 SMM-1 across the AMI band is relatively flat, suggesting 
that we seeing optically thin free-free emission. 

5.13 L1221 

The two sub-mm cores of the L1221 cloud (Young et al. 2006; 
Wu et al. 2007) were further resolved into three infra-red sources 
by Spitzer (Young et al. 2009). Centimetre radio observations of 
L1221 (Rodriguez & Reipruth 1998; Young et al. 2009) at 3.6 and 
6 cm revealed what appears to be optically thin free-free emission 



from L1221-IRS3 (SMM-1) (a = -0.16±0.54), but no detection 
in the direction of L1221-IRS1 and L1221-IRS2 (SMM-2). The 
AMI-LA 16 GHz observations, see Fig. 5 (c), clearly detect sources 
at both the location of L1221-IRS3 and towards L1221-IRS1 and 
L1221-IRS2 which are unresolved by the AMI-LA synthesized 
beam. The peak of the 16 GHz emission towards L1221 SMM-2 
is more closely associated with L1221-IRS1, which is classified 
as a Class I object, than L1221-IRS2 which is thought to be a 
YSOc. The flux density measured from the AMI-LA observations 
for L1221-IRS3 is 5i 6 = 359 ± 28 flJy, in excess of a flux density 
extrapolated from the free-free model of the centimetre measure- 
ments by a factor of almost two. The 3.3 mm measurements of Lee 
& Ho (2005) indicate that the vibrational dust spectrum of L1221- 
IRS3 does not possess a flattened tail, such as might be expected in 
the presence of a protoplanetary disc that could explain the excess 
in terms of dust emission. If, as suggested by Young et al. (2006), 
the centimetre emission observed towards LI 221 -IRS 3 is due to 
shock ionization of the molecular outflow impacting on overdensi- 
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Right Ascension (J2000) 

Figure 4. Primary beam corrected mosaic of LI 148. AMI-LA data 
are shown as greyscale and contours linearly from 3 o", where a = 
22 /iJy beam" 1 . Spitzer candidate embedded objects from Dunham et al. 
(2008) are indicated by crosses, and the positions of sub-mm cores are 
shown as unfilled squares, see text for details. The AMI-LA primary beams 
are shown as circles and the AMI-LA synthesized beam is shown as a filled 
ellipse in the lower left corner. 

ties in the surrounding envelope (Curiel et al. 1987; 1989; Anglada 
1995; Shang 2004), then it is possible that the extra flux density 
at 16 GHz is either indicative of a second shock component with a 
larger optical depth or may represent a steepening of the radio spec- 
trum due to thermal emission from the immediate envelope of the 
central source as predicted by the model of Konigl (1984). The two 
bright sources to the west of LI 221 are coincident with a known ra- 
dio source in the NVSS catalogue, the resolution of which does not 
resolve the two component sources, and we therefore identify them 
definitely as extragalactic objects, the two further sources detected 
in this field have no lower frequency radio counterpart, but also no 
known protostellar associations. 



6 DISCUSSION 

6.1 The nature of the cm-wave radio emission: the thermal 
dust contribution 

At 16 GHz we must be aware of the potential contribution to our 
measured flux densities from the vibrational, or Planckian, dust 
spectrum. Where available we use sub-mm flux densities from the 
literature at wavelengths from 350 to 1300 flm to constrain the 
thermal dust emission at 1.8 cm. In general we use flux densities 
extracted from 40" apertures, with notable exceptions being those 
objects measured by Visser et al (2002) who used a 50" aperture, 



and IRAM 04191-IRS for which 60" apertures are assumed in or- 
der to accommodate the resolution of the 1.3 mm data (Dunham 
et al. 2006). We assume optically thin dust emission with j3 = 1.5 
and fit SEDs with a single-temperature greybody spectrum: 

s?^v^ v (r tet ). (i) 

The flux densities predicted from these fits are listed in Col- 
umn [8] of Table 5, which also lists the physical parameters of the 
cores taken from the literature. The predicted values of the vibra- 
tional dust emission are in general low compared to the rms noise 
in the datasets, detectable at 3 a only in a few fields. 

6.2 The nature of the cm-wave radio emission: free-free 
emission 

Massive star formation regions have been known for a long time 
to produce radio emission, however the low luminosity protostars 
in this sample do not produce a sufficient ionizing flux to maintain 
a photoionized Hll region such as those found in hypercompact or 
ultracompact Hll. Instead it is likely that an excess of radio flux 
density relative to the thermal dust emission in these objects is due 
to either a shocked molecular outflow or a stellar wind. A com- 
pletely ionized stellar wind (Panagia & Felli 1975) was shown to 
be feasible by Curiel et al. (1989) only for objects <C 1" in size, 
as it would produce radio flux densities far in excess of those mea- 
sured, e.g. at 16 GHz for a wind of 5" angular size and assuming 
an electron temperature T e = 10 4 K, similar to Bern 48, we would 
observe a flux density of Sie ^ 50 Jy. Since the majority of our 
observed sources are extended on scales of greater than an arcsec 
we find this possibility untenable. The alternatives are that we are 
seeing radio emission from a neutral wind or molecular outflow, 
which is being shock ionized as it impacts on overdensities within 
the surrounding interstellar medium or the infalling cold dynami- 
cally collapsing medium (Curiel et al. 1987), or we are seeing ra- 
dio emission from a partially ionized and possibly collimated flow 
(Reynolds 1986). 

The cm-wave flux density from regions of shock ionization 
has been shown to be (Curiel et al. 1987; 1989) 

(Jy_\ = ,,,,,, ,,, W \( v* V' 68 

VmJyy ' VlO" 7 M yr-V Vl00kms-V 

x (kp^) (itTk) (4^) (5GHz") (2) 

Following Anglada et al. (1996; 1998) we may re-express this rela- 
tion at 16 GHz as 

/ F out \_ 2.97X10- 4 / S V D 2 \ 
Weyr^kms-V n qK j VmJykpcV' W 

where r\ = (Cl/4n) is the fraction of the stellar wind being shocked, 
F ut is the outflow force or momentum flux, equivalent to the rate 
of outflow momentum and often computed as F out = (P/ Td yn ). We 
also introduce ^(t) = t/(1— e~ T ) to allow for the fact that the 
radio emission may not be optically thin (Anglada et al. 1998), al- 
though this dependence on the optical depth, T, is very weak. 

The parameter r\ may be considered as the efficiency of the 
shock and, following Anglada (1995), a value of r\ = 1 allows us 
estimate the minimum outflow force required to explain the cm- 
wave emission as shock ionization. In Fig. 6 (a) we plot the corre- 
lation of the radio luminosity as measured at 1 .8 cm from this paper 
against the outflow force for each object where a value is available 
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Table 5. Physical properties of the AMI LA Protostar sample. 
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CB188 SMM-1 
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1.94 ±0.2, 1.5 


1,5 


3.5 C 


5 


18 


-1.795 f 


I 


L673-7-IRS 


300±100 


0.017 


0.09±0.03 


10 






21 


-2.290* 


O.VeLLO 


L723-IRS 


300±100 


0.048 












-1.703 


O.VeLLO 


L723 SMM-1 


300±100 




2,3 


1 


37 


14 


69 


-1.253 





LI 152 SMM-1 


325 




3.3 


1 


0.47 


14 


15 


-1.809* 


on 


L1148-IRS1 


325±25 


0.081 


0.23 ±0.17,0. 15 


10 








-2.120 


O.VeLLO 


L1148-IRS2 


325±25 


0.003 












-1.994* 


O.VeLLO 


Bern 48 


200 




11.05±0.20,13.5 


6 






39 


-1.609 


I/II 


L1014-IRS 


250±50 


0.087 


0.3,0.34±0.11 


10,11 


0.004 c 


7 


6 


-1.729 


O.VeLLO 


L1165-IRS 


300±50 


0.030 












-2.032 


O.VeLLO 


LI 165 SMM-1 


300±50 




11.7 


4 


1.1 


4 


33 


-1.663 


I 


L1221-IRS1 


250±50 


1.12 


2.5 


10 






32 


-1.649 


I 


L1221-IRS3 


250±50 


0.36 


0.93 ±0.16 


10 


10.1 c 


13 


29 


-1.649 






* 3 a upper limit. 

' lower limit, see text for details. 

c Value of F out from the literature corrected for D. 

(1) Furuya et al. 2003; (2) Dunham et al. 2006; (3) Andre et al. 1999; (4) Visser et al. 2002; (5) Yun & Clemens 1994; (6) Connelley et al. 2008; (7) Bourke 
et al. 2005; (8) Hayashi et al. 1994; (9) Ohashi et al. 1996; (10) Dunham et al. 2008; (11) Young et al. 2004; (12) Myers et al. 1988; (13) Umemoto et al. 1991; 
(14) Bontemps et al. (1996); (15) Armstrong & Winnewisser (1989). 



in the literature. For values and references see Table 5. We also plot 
the predicted outflow force from the Curiel model for an efficiency 
of rj = 1 (minimum required force; lower line) and for an efficiency 
of 77 =0.1. 

There are considerable uncertainties in outflow force measure- 
ments and variations in how the force is estimated. Some estimates 
are derived from the total momentum observable in the outflow 
lobes, divided by a measure of the outflow dynamical time (e.g. 
B35A SMM-1 or CB188). Others use a measure of the force de- 
rived within an annulus around the outflow origin (Bontemps et al. 
1996; e.g. L723 SMM-1, LI 165 SMM-1), which provides a better 
estimate of the force for outflows which have been only partially 
mapped. Even analyses which use similar methods can have vari- 
ations, e.g. whether the maximum velocity in an outflow is used 
to compute the dynamical time or some average/typical velocity. 
Furthermore, missing flux density (in interferometric studies), high 
optical depths of CO emission in the outflow line wings and incli- 
nation effects which reduce the observed outflow velocities all tend 
to reduce the momentum actually observed. 

Many studies follow the lead of Bontemps et al. (1996) and 
increase their outflow forces by a factor of 10 to account for op- 
tical depth and inclination effects (e.g. Visser et al. 2002). Other 
authors are able to correct for optical depth using observations of 
the rarer isotopologues of CO (e.g. L1221-IRS3; Umemoto et al. 
1993) or note that the inclination corrections will be minimal (e.g. 
L1014-IRS; Bourke et al. 2005). A consistent approach for the out- 
flows in our sample is not possible. Instead we use the values as 
stated in the literature, corrected only if we assume a different dis- 
tance (F out 0= D). For a large enough sample of protostars, the trends 
should still be apparent, albeit with a large scatter, but the force for 
any particular outflow could be under-estimated by a factor of fa 10. 

From Fig. 6 (b) we can see that three objects lie below the 



calculated minimum outflow force: L1014-IRS, L673 SMM-1 and 
B35A SMM-1. Missing flux and optical depth effects are not taken 
into account in the outflow force of L1014-IRS (Bourke et al. 
2005). However, even if the value were scaled upwards by an or- 
der of magnitude it would still fall below the predicted minimum 
force limit. A similar upwards scaling for B35A SMM-1 could 
be required, as optical depth and inclination were not accounted 
for, unlike L673 SMM-1, which has been scaled upwards follow- 
ing the prescription of Bontemps et al. (1996). If all the objects 
were treated in the same way it is likely that only L1014-IRS and 
L673 SMM-1 would remain below the minimum line. Shirley et al. 
(2007) found the same situation for L1014-IRS at 3.6 and 6 cm and 
concluded that L1014 possessed an additional ionization mecha- 
nism aside from shock ionization. This additional component is 
most likely to be provided by a self-shocked spherical wind (Pana- 
gia & Felli 1975; Wright & Barlow 1975; Reynolds 1986). A caveat 
to this explanation is that spherical winds are usually associated 
with more evolved Class II objects. 

As discussed earlier, a fully ionized spherical wind would be 
likely to produce a flux density far in excess of that which we mea- 
sure for L1014-IRS, particularly as the radio emission is slightly 
extended. An alternative is to consider the model of Reynolds 
(1986), which allows for partially ionized, collimated winds. An 
advantage of this model above that of a spherical wind is that it al- 
lows for values of the radio spectral index from an unresolved out- 
flow of a > 0.6, as has been observed in a number of cases. This 
model assumes power-law dependencies of physical parameters on 
the radial distance, r, from the central source. The first of these 
is the jet half-width, w, which is assumed to obey w = wo(r/ro) £ 
where ro represents the "ionization radius" at which the thermal 
emission begins. For well collimated flows, ro and wo are related 
by the characteristic function 6*0 = 2wo/ro 0.5. Reparameterizing 
the equations of Reynolds (1986) in terms of the physical parame- 
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Figure 6. Distribution of radio luminosity vs. outflow force, see Table 5 for 
references, (a) Data points are radio luminosities derived from the AMI-LA 
measurements at 16 GHz, see Table 5. The dashed lines show the theoretical 
relationship between 1.8 cm radio luminosity and outflow force (Curiel et al. 
1989) for an efficiency of rj = 1 (minimum required force; lower line) and 
an efficiency of 7] =0.1, see text for details, (b) Filled data points are radio 
luminosities derived from the AMI-LA measurements at 16 GHz as before; 
unfilled data points are from Anglada et al. (1998) extrapolated to 16 GHz. 
The dotted line shows a linear regression to the AMI-LA data alone, and 
the solid line shows the regression to the combined data sets. 



ters observed here we find that the predicted radio emission can be 
found using, 



5 v (mJy) = 3.2 x ft ( 



''() 



10 15 cm/ V 10GHz/ 1 V 10 GHz 



(4) 



where 



10 GHz 



3.7 x 10" 'ft 



lO-SMeyr-'kms" 1 



ro \"3/ V, VY T, 



10 15 c m y V200kms - V \10 4 K 



(2.1) 2 



q T (a-2){a + 0.\y 



(5) 



(6) 



and 



c=(1 + £ + 9t)/?t : 



a - 



2.1 



(7) 



The parameters qj and q T , corresponding to the notation of 
Reynolds (1986), represent the power-law indices for temperature 
and optical depth, respectively. The values of these parameters for 
various models may be found in Table 1 of Reynolds (1986), as 
may the predicted value of a. The sim dependence in this formula- 
tion takes into account the inclination angle, i, of the outflow to the 
observer's line of sight, xq is the fractional ionization of the gas, x, 
at the ionization radius, which is a tunable parameter in the model 
and affects the optical depth through the free-free absorption coef- 
ficient , kv, in the standard way: T °c x 2 . A caveat to this model is 
that there is no explanation of how the ionization is produced, and 
a value must simply be assumed. 

By combining Equations 4 and 5, we can see that this model 
predicts a correlation between the radio luminosity and the outflow 
force of, 



. 2.1/2(2-a) 
Had 



(8) 



For instance, in the case where a = 0.6 this correlation becomes 

From the number of data points in the AMI-LA data alone 
which have measurements of F out , see Fig. 6(a), there is no evi- 
dence of a correlation and the Pearson product momentum corre- 



lation coefficient is low: r. 



XV 



0.12. To enhance the trend we also 



include datapoints extrapolated from the measurements of Anglada 
et al. (1998), Fig. 6 (b). These data are extrapolated from measure- 
ments at 3.6 cm using a power law with spectral indices derived 
from data at 3.6 and 6 cm. The flux density measurements, derived 
spectral indices and measured outflow forces for these data may be 
found in Table 5 of Anglada et al. (1998). 

Including the additional data, the correlation between outflow 
force and radio luminosity is still weak (r xy = 0.40). This may be 
due to a lack of data, but could also be a consequence of the uncer- 
tainties in the quoted values of F ou t. The correlation we find from 
the AMI-LA data alone is 



log [F OTt (M yr 'kms 



= -(4.48 ±1.26) 

+(0.26±().74)log[L 1 . 8cm (mJykpc 2 ) 



(9) 



If we also include the extrapolated data from Anglada et al. (1998) 
we find 



log[iw(M yr 'kms ')] 



-(3.54±0.58) 

-(0.73 ± 0.38) log[Li. 8cm (mJykpc 2 ) 



(10) 



see Fig. 6(b). This correlation is very similar to that found by 
Shirley et al. (2007) at 3.6 cm. 

Due to the short frequency coverage of the AMI-LA band and 
the low signal to noise of the objects in this sample we consider the 
errors on any spectral index derived from AMI-LA data alone too 
large to allow a reliable quantitative analysis. However, six of the 
objects in our sample have additional radio measurements at lower 
frequencies, these are listed in Table 6, and we use these to exam- 
ine the spectral behaviour of the radio emission from these cores. 
Where the source has not been detected we take a 3 a upper limit 
on its flux density, with a equal to the rms noise as quoted in the 
reference. We find that all of the sources have rising spectral in- 
dices, indicating that they all possess partially optically thick free- 
free emission, with the exception of L723 SMM-1. This source has 
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Table 6. Radio Spectral Indices. 



Name 


S(6cm) 


5(3.6cm) 


5(2(1.8) cm) 


a 


refs. 


IRAM 04191 


480 ±30 


140 ±20 


160 ±50, 172.9 ±17 


0.45 ±0.20 


1, this work 


L1521F 




<66 


141 ±16 


> 1.10 


2, this work 


B35A 


< 120 




943 ± 67 


> 1.77 


3, this work 


L723-SMM1 




740 ±42 


690 ±29 


-0.10±0.21 


4, this work 


L1014 


88±11 


111 ±8 


299 ±21 


0.71 ±0.11 


5, this work 


L1221-IRS3 


192 ± 27 


177 ±24 


359 ±19 


0.83 ±0.10 


6, this work 



[1] Andre et al. (1999); [2] Harvey et al. (2002); [3] Rodriguez et al. (1989); [4] Anglada et al. (1991); [5] Shirley et al. (2007); [6] Young et al. (2009). 



a spectral index consistent with optically thin free-free emission. 
We also note that the outflow force quoted for L723 SMM- 1 is an 
order of magnitude larger than that of any other source (although 
this incorporates the upwards scaling of Bontemps et al. (1996)). 
Such behaviour may be consistent with the shock ionization taking 
place further from the central protostar, since the model of Curiel 
et al. (1987; 1989) predicts that the distance of the shock from the 
star, R, is proportional to the outflow force, 

fAUi.22xl0 3 T-°- 5 f ^? OD 

\A\]J VMayr-ikms-V 

(Anglada et al. 1998). However, the more detailed treatment of 
Reynolds (1986) shows that although the optical depth is a power 
law function of the distance from the central star, t(r) = To(r / Vo) 9t , 
the index of the power-law can have a wide range of values, 
— 4 < q x < —0.9, depending on the morphology of the outflow and 
whether recombination takes place. Substituting these values into 
Equ. 1 1 we find 

(R_\ ( F out y*»<*n 
\m) ^Meyr-ikms- 1 ; 

Given that this is an oversimplification, it is therefore perhaps un- 
surprising that we find no correlation between spectral index and 
outflow force (r = 0.12). In the case of the standard spherical flow 
(Panagia & Felli 1975; Wright & Barlow 1975) q t = -3, which 
would produce a correlation of R °c F ou l * might then interpret 
the observed lack of correlation, both in this work and as found 
by Shirley et al. (2007), as evidence that the majority of cm-wave 
emission from protostars is not produced, at least solely, via this 
model. 



6.3 The Very Low Luminosity Radio Correlation 

The correlation between radio and bolometric luminosity for lower 
mass protostars was first described for 3.6 cm emission by Anglada 
(1995) and recently updated by Shirley et al. (2007) who also in- 
cluded a correlation for 6 cm emission. This correlation is thought 
to arise as a consequence of the correlation between outflow force 
and bolometric luminosity (Shirley et al. 2007) as the degree of 
shock ionization is expected to be enhanced for more powerful out- 
flows from higher luminosity sources (Anglada 1995). From our 
1.8 cm data we find correlations of 

loglLLgc^mJykpc 2 )] = -(1.81 ±0.29) 

+(0.54 ± 0.38) log[L bol (L )], (13) 

log[L L8cm (mJykpc 2 )] = -(1.23±0.65) 

±(0.59 ±0.57) log [L IR (L G )], (14) 



with correlation coefficients of r xy = 0.77 and 0.84, respectively. 
These regressions are fitted excluding data that are limits only. 
The correlation of bolometric luminosity with 1.8 cm luminosity 
is similar to that found for 3.6 cm data by Shirley et al (2007) 
who found log[L 3 . 6cm (mJykpc 2 )] = -(2.24 ± 0.03) + (0.71 ± 
0.02) log[L|, o i(L0)]. We note that there appears to be a trend in both 
the normalization and multiplicative factors in the correlations as a 
function of wavelength, although we note that the limited number 
of wavelengths sampled may make this trend misleading. 

The correlation between radio luminosity and Li R provides a 
way of investigating the correlation with the internal luminosity of 
these sources. Dunham et al. (2008) define VeLLO sources as those 
possessing an internal luminosity L; nt < 0.1 L and use an approx- 
imately constant ratio of L mt /Lyg_ = 1.7 to distinguish objects with 
^ir ^ 0.05 L0 as VeLLOs. From Fig. 7 (a) we can see that 8 of our 
sample satisfy this criterion. Since L; nt is not a directly observable 
quantity and may only be found through modelling the source, we 
can use the radio luminosity as a proxy for the internal luminosity 
of these objects by re-expressing the correlation in Equ. 14 as 

L int (L ) ~ 112x [LL^mJykpc 2 )] 1 - 7 , (15) 

or, following Dunham et al. (2008) and normalizing to a distance of 
200 pc, we can further define the internal luminosity radio estimator 
as 

L£?(L S ) = 1.3 x lO-^.gcnGuJy)] 1 - 7 . (16) 

We do not consider this estimator to be as accurate as the 
70 Jim flux estimator (Dunham et al. 2008) as it relies on both the 
relationship of L lnt to L [R which is currently modelled as a simple 
ratio, and also because the correlation of L\ % cm with L IR is not as 
tightly defined as that of F70 to L mt . Two pieces of further work 
are required to improve the radio estimator of internal luminosity 
(1) a better constrained relationship between L; nt and L[ R , which 
uses additional data to investigate and fit the low luminosity ex- 
cess of Li nt /L IR above the constant ratio; and (2) more radio data 
to tighten the observed correlation between Lj R and L ra( j. However, 
it does provide a simple method to estimate the approximate evolu- 
tionary state of an object in the absence of complex models. 



6.4 The Non-Detections 

There are four protostellar (i.e. not starless) objects in our sam- 
ple which are undetected by the AMI-LA at 1.8 cm. The first of 
these is IRAS 04191. From our measured correlations both the L IR 
and L(, i of this object should place it firmly above the detection 
threshold in our map. From the correlation of radio and bolometric 
luminosity we would predict a flux density of S\% cm ~ 518/iJy. As 
IRAS 04191 lies approximately one- third of the way to the primary 
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beam FWHM in our image we would expect to see IRAS 04191 de- 
tected at ~ 25 a in Fig. 1 (a). However, IRAS 04191 is an evolved 
Class I object and if we make the assumption that the cm-wave 
emission we are seeing is a consequence of shock ionization due 
to a molecular outflow, we may explain the lack of radio emission 
from IRAS 04191 by the fact that Class I objects can have outflows 
from an order of magnitude to several orders of magnitude weaker 
than those of Class O sources (Bontemps et al. 1996; Curtis et al. 
2010). The lack of cm-wave emission towards IRAS 04191 is in 
this case perhaps not surprising and may simply be due to a very 
weak outflow. 

The second of our non-detections is L673-7-IRS. This object 
was identified as an embedded protostar by Dunham et al. (2008) 
from the Spitzer c2d data and has L|, i = 0.09 ± 0.03. It was re- 
cently confirmed as a protostar through deep CO(2— >1) observa- 
tions (Dunham et al. 2010) which measured an outflow momen- 
tum flux of F out > 2.1 x 10~ 6 M e yr" 1 kms" 1 (or F out > 4.5 x 
10~ 6 M yr~' km s _1 using the force per beam method (Fuller & 
Ladd; Hatchell et al. 2007). We can see from Fig. 6(a) that with an 
outflow of such low momentum flux any associated radio emission 
may lie below our detection limit, possibly by more than an order 
of magnitude. 

L1152-IRS is the third of our non-detections. The AMI-LA 
map of L1152-IRS shows some small structure at the pointing cen- 
tre, coincident with the position of the sub-mm core. It is likely that 
the increased noise level in this dataset due to satellite interference 
has prevented us from making a significant detection of this object. 

The final non-detection in our sample is LI 148-IRS2, and in- 
deed L1148-IRS1 may be considered a non-detection also. L1148- 
IRS2 is the very lowest LrR object in our sample and was classified 
by Dunham et al. (2008) as a Group 6 source under their classifi- 
cation scheme. Group 6 has the lowest probability of being a true 
protostar and it is possible that L1148-IRS2 is in fact an external 
galaxy. However, for a true source of such low predicted Lir as 
L1148-IRS2, we can see from the radio-IR-luminosity correlation, 
Fig. 7, that the AMI-LA observation would need to be at least an 
order of magnitude deeper in order to detect this object at 3 <J. In 
the case of L1148-IRS1 it is possible that the Spitzer spectrum is 
confused by an alternative physical process. It has been suggested 
(Hatchell & Dunham 2009) that the Spitzer spectra of some low 
luminosity objects may confuse true embedded sources with H2 
shocks or PDR excitations (Smith et al. 2006). In this case such a 
hypothesis is supported by the lack of a compact sub-mm core in 



Table 7. Detection Statistics. 



Class 


Total 


Det. 


Non-det. 


Rate 

(%) 


VeLLO 


8 


5(2) 


3(3) 


63(100) 





4* 


3 


1* 


75 


I 


8* 


6 


2* 


71 


11 


1* 


1* 




100 


Starless 


3 





3 






the SHARC-II 350 flm observations of this field, which place an 
upper limit on the sub-mm flux density of S^so^m < 0. 1 1 Jy (Wu 
et al. 2007). Although a compact peak in the sub-mm does not con- 
clusively indicate the presence of protostellar activity, it is difficult 
to reconcile such activity in the absence of one. 

6.5 The Detections 

The detections and non-detections are summarized in Table 7. The 
bracketed values indicate the number of VeLLOs within this bin 
that have not been confirmed independently of their Spitzer identifi- 
cation and which may be confused by the physical mechanisms dis- 
cussed in § 6.4. A starred value indicates a number which includes 
an object with a borderline class (e.g. LI 152 SMM and Bern 48). 
The detection rate does not vary a great deal across the protostellar 
objects. A detection rate of zero towards starless objects confirms 
that the radio emission we see arises from a physical process other 
than a flattened dust tail due to vibrational emission. This highlights 
the ability of 16 GHz radio observations to identify protostellar ac- 
tivity in cores which do not possess definitive sub-mm or infra-red 
measurements. 



7 CONCLUSIONS 

The rising radio spectrum of the cm-wave emission from protostel- 
lar objects makes higher radio frequencies more suitable for their 
detection than the longer wavelengths often employed. The low in- 
tensity of thermal dust emission at frequencies below ~ 50 GHz 
greatly reduces the possibility of confusing the observed radio 
emission associated with molecular outflows or stellar winds, and 
therefore an indication of protostellar activity, with the cold dust 
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tail from starless cores. The combination of these two effects makes 
16 GHz an effective frequency at which to distinguish protostellar 
objects and starless cores. This situation also extends to the very 
low luminosity objects, VeLLOs, and embedded objects since the 
dense cores surrounding these objects are transparent to the longer 
wavelengths of the radio emission. 

Two of the objects detected here, L723-IRS and L1165-IRS, 
were classified as poor candidates for embedded objects by Dun- 
ham et al. (2008). However the cm-wave emission we detect to- 
wards these sources might suggest that some protostellar activity is 
present. As previously noted these objects may be external galaxies, 
or possibly later type stars which have been mis-identified in the 
MIR data, both of which could have radio counterparts. However, 
in both cases the morphology of the radio emission is not point-like, 
notably so in the case of LI 165-IRS, which suggests that the radio 
emission is most probably not from a stellar wind in a later type 
star as this is generally well modelled by a spherical morphology 
and would appear point-like to the AMI-LA synthesized beam. If 
the radio emission arose from an external galaxy we would expect 
to see a steep spectrum synchrotron component and therefore also 
be able to identify the source in lower frequency radio data. More 
information is necessary to conclusively identify these objects as 
protostellar cores. 

The radio emission from the sources detected here is gener- 
ally <C 1 mjy and has a relatively compact morphology, although 
CB188 SMM-1 is possibly an exception to this situation. As such 
it should not provide a major source of confusion for spinning dust 
studies of the surrounding clouds that contain these cores, which 
are generally conducted on much larger angular scales. 

The results presented here have shown that, in the absence of 
70 jllm data, it is possible that the radio luminosity may be used as a 
proxy for the internal luminosity of embedded protostellar objects 
and to enable a larger sample of VeLLOs to be identified. How- 
ever, at present the correlation presented here suffers from a lack of 
data and will require further observations to define the trend more 
completely. 
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